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ABSTRACT 


Changes in ambient temperature and humidity affect the 
exhaust emissions of a gas turbine engine. The .results of 
a test program employing a JT8D-17 combustor are presented 
which quantize the effects of these changes on carbon 
monoxide emissions at simulated idle operating conditions . 
Analytical results generated by a kinetic model of the 
combustion process , and -reflecting changing ambient 
conditions , and a comparison of these with the experimental 
results are also given. It is shown that for a complete 
range of possible ambient variations, significant changes 
do occur in the amount of carbon monoxide emitted by a 
gas turbine at idle, and that the analytical model is 
reasonably successful in predicting these changes. 
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. INTRODUCTION 


Aircraft gas turbine engines utilize hydrocarbon 
fuels, which are allowed to burn inside the combustion 
chambers of the engine, resulting in heat release-, which 
is converted to motive power by the turbine and the- nozzle. 
The bumiing process also results in the predominant 
products of hydrocarbon-air combustion, namely, carbon 
dioxide and water. However, the combustion process is 
not ideal, and small amounts of other products are also 
formed. These auxiliary products may be -categorized as 
the products of incomplete combustion,. oxides of nitrogen, 
and -products arising out of .impurities in the fuel. 

The products of incomplete combustion include carbon 
monoxide, unburnt hydrocarbons and smoke; the oxides of 
nitrogen include nitric oxide and nitrogen- "dioxide , 
collectively called NO ; impurities in the -fuel give rise* 
to o-ther products .such as oxides of sulfur. 

Smoke and NO emissions are found to occur in quantity 
when the engines are operating under high power settings . 
Smoke has been identified as being produced in the fuel*^ 
rich regions of the combustor, and by proper combustor 
design, taking care to avoid such regions, smoke emission 
has been significantly reduced. Production of oxides of 
of nitrogen in quantity arises out of extended contact 
between oxygen and nitrogen in the air at high temperatures 
within the combustor. 
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Carbon monoxid© and unburnt hydrocarbons are pro~ 
duced in quantity under low power settings of the engine. 

At such settings , both the mass flows of air and fuel are 
low, but the mass flow of air is much higher than that of 
the fuel, resulting in very lean primary zones. Carbon 
monoxide and unburnt hydrocarbons result from such zones, 
because the flame is reiatively cold at the low fuel/air 
ratios encountered here, due to which, oxidation of 
carbon monoxide and hydrocarbons is inhibited. 

These emissions, namely, oxides of nitrogen, smohe, 
carbon monoxide, and unb.urnt hydrocarbons have been -called 
"pollutants” by the United States Environmental Protection 
Agency. While -the actual amounts of these pollutants 
discharged from a single engine may be diminutive, it 
is important to note that these emissions occur entirely in 
the vicinity of air transportation centers , at the ground 
level, and that thousands of engines may be operated in 
the course of a single day. Engines operate under low 
power settings during idling and taxiing , and under high 
power settings during take off. The consequent release, 
.into the local environment, of large amounts of these 
pollutants is of great conceam. 

Programs are in progress to assess the actual impact 
of such emissions from aircraft gas turbine engines, and 
other sources, on the earth environment, on the climate 
patterns, on the inhabitants, and on vegetation. Xt is 
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beyond do^±)t that there are both localized and dispersed 
effects . In addition to just vitiating the atmosphere, 
there is the danger of setting in motion other reactionary 
processes which may result in irreversible consequences. 

The global impact of pollution on the climate is of great 
concern, but is yet to be recognized. However, short-term 
and localized effects such as smog, effects on vegetation 
and -adverse physiological effects on humans are well-known. 
The danger of depleting the ear-th's protective ozone layer 
in .the stratosphere .and consequently raising the levels of 
undesirable ultra-violet radiation at the earth's surface, 
cuasing physiological problems such as skin cancer, is also 
of great concern. This is an. example of the possible long- 
term effects of indiscriminate pollution, and indicative 
of the gravity of the situation, if the pollution problem 
is not taken well in hand, at .an early stage. 

Towards this goal, namely, to limit pollution from 
future powerplants, the USEPA made a detailed study of the 
situation and issued a set of regulations. For aircraft gas 
turbine engines, in particular, standards "were issued in' 
1973^^^. These standards were in the form of EPA parameters 
or EPAP's, and were given as the emission -index over a 
landing/take-off cycle. The emission' index is defined as 
the grams of pollutant released per kilogram of fuel 
consumed. Engines , whose pollution levels exceeded these 
limits would not be certified for operation. Also, the 
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EPAP ' s . related to pollution levels -at standard day 
conditions of 1 atmosphere pressure., 289 “K temperature, 
and zero h\amidity (the standard humidity has been changed 
to 60% relative humidity) * Since actual atmospheric 
conditions rarely meet the standard day criteria, and since 
emission index measurements made for certification purposes 
may vary, depending on the particular ambient conditions 
of measurement, it is necessary to. have a method. for 
correcting non-standard emission index measurements for 
ambient condition variations . Towards this end a correction 
factor was defined, as follows: 


Correction Factor 


Emission Index under Standard Conditions 
Emission Index linder Nonstandard Conditions 


Knowing the nonstandard conditions of measurement and thus 
the relevant correction factor, a priori, and the measured 
emission index (nonstandard) , simple multiplication of 
the two. will result in the emission index which would have 
occurred (for the particular engine) , had the conditions of 
measurement been standard atmospheric day conditions . 

Additionally, it is desirable to have a general equation, 
at least applicable to one particular type of engine, for 
the correction factor, incorporating ambient and other 
effects. This and similarity between different ' engine 
types would greatly ease the regulatory task of enforcing 
the EPA standards. Correction factor correlations, including 
ambient temperature^ pressure , and humidity effects are 
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therefore desired. 


Quite early, Lipf ert ' ^ noted the effects of changes 

in inlet pressure, temperature, and hmidity on the emission 

levels of oxides of nitrogen, produced by an engine 

operating under take-off power conditions. Si^bsequently , 

numerous correction factor correlations were formulated. 

A compilation and evaluation of these and other correction 

( 4 ) 

factors have been given by R-ubins and Marchionna 

. ( 5 ) 

For a combustor operating at idle, Marzeski and Blazowski 
developed additional correlations for the effects of 
nonstandard ambient pressure and temperature on all emissions. 
For production samples of a given engine, the effect of 
ambient temperature and pressure on all emissions over 
the complete thrust range has been correlated by Sarli 
et al. . 


With the exception of some limited engine test results 

f 7 ^ f o > 

given by Nelson^ ^ , Hosier et al. , and the work reported 

(9 \ 

by Allen and Slusher , the effect of humxdity on carbon 
monoxide emissions from idling gas turbine combustors has 
received little attention, although the extreme sensitivity 
of carbon monoxide oxidation to the presence of water 
vapor is well-known. 

To ascertain the effects of ambient conditions , humidity 


in particular, on carbon monoxide emissions from idling 
gas turbine combustors , a two-pronged research effort was 
initiated, encompassing both experimental and analytical 
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work. Experimentally, a combustor rig was employed to 
simulate changing combustor inlet conditions as generated 
by changing ambient conditions and carbon monoxide 
emission levels were measured at the combustor exit for 
each simulated ambient condition. The results revealed 
that ambient hiomidity and temperature can, indeed, have 
significant effects on carbon monoxide emission levels. 

As far as the analytical effort is concerned, it is 

well known that equilibrium considerations lead to 

erroneous results , in that the predicted carbon monoxide 

levels can differ from measured levels by several orders 

of magnitude. Hence, chemical kinetic considerations may 

be more appropriate in describing the combustion process 

in a gas turbine combustor. It is also reasonable to 

construct a combustor model, with provisions for incorporating 

both local combustor conditions and ambient conditions , in 

order to be able to compute the carbon monoxide -emission 

levels with respect to changes in ambient conditions . Some 

details relevant to carbon monoxide production in gas 

turbine combustors have been given by Morr et al. , and 

a less detailed model, but including limited ambient 

effects has been presented ny Sarli et al. ^ . A review 

of various gas turbine combustor modelling concepts .and 

( 11 ) 

emissions has been given by Osgerby , wherein sophrstxcated 
combustor models employing perfectly stirred reactors , 
partially stirred reactors and combinations of these and 
plug-flow reactors, such as in what is known as a 


6 



" Swithenbank’ model” are also discussed. 

In attempting to define a tractable, yet accurate 
model of the combustor, conflicts will arise. It must be 
borne in mind that even sophisticated combustor models 
have not been totally satisfactory in predicting emission 
•levels. I!,urther, the increased complexity and sophisti- 
cation in computation techniques , necessary .when such 
models are employed, may, in fact, obscure the principal 
features of pollutant production, and, also result in high 
computation costs . 

Plug-flow reactor models, in comparison, are known 
for their simplicity and ease of computation. Reasonable 
results have been obtained using plug-flow reactors in 
combustor modelling ^ (1^) ^ 

of simplicity, ease of computation, reduced cost of 
computation, and retention of the principal features of 
pollutant production, plug-flow reactor modelling is 
employed in this research effort. In the analytical 
program, a chemical reaction kinetic scheme was applied in 
the combustor modelling, where initial species concentra- 
tions reflected both combustor characteristics and changin< 
ambient conditions , 

The goals of the present work are to present .a -simple 
combustor model and to apply this model to predict the 
effects of ambient temperature and humidity on the carbon 
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monoxide emission levels from an idling gas turbine 
combustor, to develop correction factors for these effects 
and to compare the calculated and measured emission 
levels and correction factors. 
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EXPERIMENTAL PROGRAM 


The aims of the experimental program are to measure 
carbon monoxide emission levels from gas turbine combustors 
operating under idling conditions for various simulated 
ambient temperatures and humidities, and to ascertain 
their effects on the emission levels. 

The experimental program was conducted at the NASA 
Lewis Research Center, employing their closed-duct test 
facility, which is described in detail by Fear I This 
facility consists of a combustor' rig , with' provisions for 
supplying appropriate amounts of humidified non-vitiated 
air at .conditions which simulate various compressor 
pressure ratios, at the combustor inlet. 

'A single JT8D-17 combustor can was installed in the 
duct, along with the appropriate instrumentation .for 
measuring static pressure, total -temperature and gas chemi- 
cal composition as shown in Pig. 1. The can was' supplied 
with appropriate amounts of non-vitiated air and Jet A 
fuel to maintain a constant overall fuel/air ratio. 
Demineralized water was injected well upstream of the 
combustor inlet, to insure complete vaporization, in order 
to simulate ambient h'umidity. The temperature, pressure, 
and humidity at the combustor inlet were controlled to 
reflect desired ambient conditions, pressure ratios, and 
compressor efficiency. 
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Humid air and fuel were allowed to undergo combustion 
in the JT8D-17. can under these simulated conditions, and 
the total temperature, static pressure, and chemical 
composition were measured to SAE .specifications at 
the various locations shown in Fig, 1. 

Test Conditions 

The idle operating conditions , both nominal and as 
tested, are given in Table I,. In relating ambient variables 
to combustor inlet variables, compressor pressure ratios 
of 2, 3, 4 and 5 were chosen, along with a compression 
efficiency of 80 percent. The mass flow through the. 
combustor consists of both .air and water vapor, "the com- 
bination of which may be 'regarded as the oxidizer. The 
fuel flow was controlled to maintain a constant_^f_uel/air 
ratio- and not a constant fuel/oxidizer ratio. The actual 
mass flow through the combustor was calculated on either 
a constant compressor discharge Mach number basis, or a 
constant reference velocity basis. Also, since the 
combustor geometry was fixed, three different values of 
overall fuel/air ratio were run, in order to affect local 
fuel/air ratios within the combustor. 

The combustor inlet pressure, temperature, and humidity 
were controlled to reflect desired ambient conditions of 
244, 289 and 322°K temperature and 0, 50 and 10.0 percent 
relative hiamidity, except at 244 ®K, at which temperature 
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an extremely small amount of water vapor corresponds to 
saturation, so only one value of relative humidity could 
be obtained, which was approximately 60 percent. The 
water content of the air supplied by the preheater was 
continually monitored and was quite small. 

Experimental Results 

Representative values of the measured carbon monoxide 
emissions from the JT8D-17 combustor can are given in 
Figs . 2 , 3 , 4 and 5 , in terms of the carbon monoxide 
emission index against the average combustor discharge 
temperature . 

Each separate figure corresponds to one simulated 
compressor discharge condition - a pressure ratio of 3 
or 4, and either a constant compressor discharge Mach 
number or a constant reference velocity. On' each figure, 
three sets of data are presented, one for each overall 
fuel/ air ratio of 0.015, 0.011 and 0.007. 

Identically-shaded symbols correspond to one particular 

fue.l/air ratio, differently flagged symbols depict the 

effect of ambient temperature, and identically flagged 

different symbols depict the effect of h\m.idity. The 

various pressure ratios, compressor discharge conditions, 

and appropriate symbol legend are given in the figures. 

The effects of ambient temperature and humidity are felt 

both on the combustor discharge temperature and the carbon 

( 

monoxide emission index as follows; 
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_ As' .regards • the 'combustor discharge temperature: 

(i) For all other conditions remaining unaltered, a 
larger fuel/air ratio gives a higher coxnbustor discharge 
temperature (compare differently shaded' symbols) . 

(ii) For a fixed' fuel/air ratio, and zero humidity, the 
combustor discharge temperature increases as ambient 
temperature is increased {compare identically shaded 
circular symbols) . 

(iii) -For a fixed fuel/air ratio and ambient temperature, 
an increase in the relative humidity leads to a. decrease 

in the combustor discharge temperature (compare identically 
shaded, identically ‘flagged symbols) . 

As regards the carbon monoxide emission index: 

■(i) For a fixed set of ambient conditions, an' '■•increase 
in the fuel/air ratio leads to a decrease in the carbon 
monoxide emission index (compare differently shaded symbols) . 

(ii) For a fixed fuel/air .ratio and .zero humidity, .an 
increase in the ambient temperature causes a decrease in 
the emission index ■ (compare identically shaded, differently 
flagged circular symbols)., 

(iii) For a fixed fuel/air ratio and- given ambient tem- 
perature, an increase in the relative hxamidity leads to .an 
increase in the emission index, an effect especially 
noticeable at the highest ambient temperature of 322°K, where 
saturation corresponds to 8.12 percent by mass of water 
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vapor (compare identically shaded, identically flagged, 
symbols.) . 


(iv) 


(V) 


For all fue'l/air ratios-, the slope- 


9EI 


CO 


ST. 


) 




is nearly identical. 
. 8EI, 


The slope 


'CO 


ST, 


) 


■4 FH=100% 

fuel/air ratio . 


increases with decreasing 


(vi) At the highest fuel/air ratio, Ct.015, the slopes 


.-r^) and 

-4 . RH=0% 


) are nearly identical. 

4 KH=100% 


As regards the basic trends noted above, there is 
little difference between a constant reference velocity 
condition or a constant compressor discharge Mach number 
condition: However, a comparison of figures corresponding 

to pressure ratios of 3 and 4 reveal that a decrease in 
the pressure ratio, for both the abpve conditions, leads 
to an increase in- the carbon monoxide emission 'index. 
Emissions, data collected at other simulated pressure 
ratios reflect similar trends -with regard to changing 
.ambient conditions. . 

Other combustor data is also surprisingly similar. 

. . - (5) 

The emissions data of Blasowski et al. ' ' for a T-r56 

combustor, employing two different fuels under idle con- 
ditions , but with different primary zone characteristics 
is presented in Pig. 6 -in a form similar to the JT8D-17 
data. During these tests, the geometry of the combustor 
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was modified slightly, to allow adjustment of the primary 
-zone fuel/air ratio, while keeping the overall - fuel/air 
ratio unaltered. The relative humidity was close to zero 
and -the simulated compressor pres'sure ratio was three. 
Although there are differences in the actual carbon 
monoxide emission index produced by the JT8D-17 and the 
T-56 -combustor under idle conditions, the slope: 

f 

RH=0% 

in both figures 5 and 6 , depicting the anbient temperature 
effect is nearly identical, especially for the lean primary 
zone of the T-56. Such a similarity between -different 
combustors, if discovered-, would considerably ease the 
regulatory task of developing corrections .for nonstandard 
ambient conditions , 
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ANALYTICAL PROGRAM 


The experimental results indicate "that carbon monoxide 
emission levels are decreased by an increasing fuel/air 
ratio, pressure ratio, and ambient temperature, while 
they .are increased by an increasing ambient humidity. In 
order to understand this behavior it would seem reasonable 
to construct a combustor model. Several factors are of 
importance when such a model is desired, in order to be 
able to predict pollutant emission levels from gas turbine 
combustors. A typical gas turbine combustor consists 
of a casing to contain the air flow, within -which, there 
iS' a means for mixing fuel and air, and burning, the' 
mixture -under stable conditions.. The actual combustion 
process is extremely complex, -and includes recirculation 
of burnt products, incomplete mixing, delay in mixing and 
burning, quenching by the walls of the combustor, and 
further, these processes occur simultaneously. Detailed 
investigations by Mellor et al. of gas turbine- com- 

bustors have revealed the existence of such complex 
processes in the confines of combustors. 

A realistic computational model of the combustor 
would have to incorporate all, these complexities , and in 
addition, have provisions for taking into account the 
effects of ambient conditions on the operating parameters , 
such as pressure, temperature, and chemical composition. 
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However, it has not been possible to include these effects 
and have a reasonable computation time, and hence simpler 
models are highly desirable. 

The combustion process in a gas turbine combustor is 
known to generally occur in three principal zones.: the 

primary, secondary, and dilution. These zones are shown 
for a typical gas turbine combustor in Fig. 7. The primary 
zone consists of a fuel injection region and a highly 
stirred region. Most of the combustion, to the extent 
of 90 percent, occurs in this zone- The seconda3:y zone 
completes the coinbustion to about 99 percent and initiates 
dilution, which is completed in the dilution zone, to 
provide an acceptable turbine inlet temperature. All 
combustor models so far suggested invariably -subdivide 
the combustor into these three zones , and differ basically 
in the manner in which the primary zone is 'modeled, 
consideration of droplet evaporation and combustion, and in 
the chemistry of pollutant production. A detailed dis- 
cussion of these aspects is given in a review paper by 
( 11 ) 

Osgerby 

In the model presented in this research effort, it 
is suggested that each such zone within the combustor be 
treated as a plug-flow reactor. In a plug-flow reactor 
of the type employed in this analysis, there is homogeneous 
chemical reaction between perfectly mixed fuel and 


16 



oxidizer, under conditions of constant pressure and 
temperature. The pressure here depends on the ambient 
pressure and compressor pressure ratio, and the tempera- 
ture is taken as the adiabatic flame temperature for each 
particular ambient condition and local fuel/air ratio. 

In actuality, adiabatic flame temperatures are not 
attained in the combustion zones due to the fact that 
chemical equilibrium is not reached in the combustor. 

The chemistry in such plug flow zones may be described 
adequately by a set of chemical reactions, or kinetics. 

Since the actual kinetics representing the oxidation of a 
hydrocarbon fuel such as Jet A are only poorly understood, 
a simpler fuel, methane, whose oxidation kinetics are 
well documented, is employed. 

It is in keeping with one of the goals of this research 
effort, namely, simplicity, that plug-flow reactors are 
being used to model the conibustor. Plug-flow models ignore 
all the complexities in the flow previously mentioned, 
and provide an idealization of the combustion process . 
However, they need not be considered as over-simplifying ' 
the matter, because even relatively 'sophisticated combustor 
models have not been totally satisfactory in accounting for 
pollutant emission levels. Some of these sophisticated 
models are presented and discussed in detail by Morr et 
al.^^®^, Osgerby^^^^, Sheppard^^^\ Pratt et and 
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( 21 ) 

Boccio et al. . Plug-flow models on the other hand, 
have been used with limited success by Seery and Bowman 
Marteney , and Starkman et al. Further, complex 

combustor models require sophisticated computation 
techniques which indeed obscure the principal features 
of pollutant production. In the interests of simplicity 
and reduced computation time and costs, plug flow 
modelling is employed in this research effort. 

As mentioned earlier, plug flow reactors allow 
perfectly mixed fuel and oxidizer to undergo chemical 
reaction in a constant pressure, constant temperature 
environment . The temperature • in a plug flow reactor 
corresponds to the adiabatic flame temperature which is 
calculated using the NASA CEC- 71- computer program ^ ^ . 

The computer program requires, as input, the reactant 
composition, pressure, and initial temperature. The 
reactant composition reflects the ambient humidity, and 
the initial temperature reflects the ambient temperature, 
the pressure ratio, and the compressor efficiency. The 
reactant inlet temperature, (i.e., at combustor inlet) 
is given by : 

Y-1 

T- = T [1 + ^ {-(PR) - 1}] (1) 

o o ’‘c 

For an explanation of the symbols, please refer to the 
" N ome nc 1 a ture ” . 
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The initial chemical composition of the 

reactants , 

which consist of fuel 

. and oxidizer, are 

calculated as 

follows : 




Oxidizer: (Proportions by mass) 



Nitrogen 

(N^) 0.7555 



Oxygen 

(O 2 ) 0.-2313 - 



Argon 

(Ar) 0.0127 



Carbon Dioxide 

(COo) 0.0005 


(2) 

Dry Air 

1.0000 



Water Vapor (H 2 O) = 

Dry Air x W^ x RH 




1 X W X RH = W 

X RH 

(3) 


s s 



Oxidizer total mass 

= Dry Air + Water 

Vapor 



== 1 + (W X RH) 

•9 


(4) 

The mass fraction of 

each component in the oxidizer is 

given by : 




Nitrogen 

= 0.7555 /(I + W 

s 

X RH) 

/ 

Oxygen 

= 0.2313 /(I + W 

s 

X RH) 

r 

Argon 

= 0.0127 /j:i +. w . 

X RH) 

f 

Carbon Dioxide 

= 0.0005 /{I + W 

X RH) 

, and 


s 


(5) 

Water 

= W X RH/(1 + W 
s s 

X RH) 

• 
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The fuel consists of methane which is given by: 


Methane (CH^) 


Dry Air x fuel/air 
1 X f/a - f/a 


( 6 ) 


Also, Methane = Wet Air x fuel/oxidizer 

= (1 + W^ X RH) X (f/o) (7) 

The relative amounts by weight, of the various components 
comprising the reactants and their initial states can now 
be calculated, if the fuel/air ratio or fuel/oxidizer ratio 
and the ambient conditions and pressure ratio are specified. 
For each such specific condition, the NASA CEC-71 computer 
program calculates the adiabatic flame temperature • and 
equilibrium chemical composition. 

The kinetic chemical calculations in the plug flow 
reactors is carried out using the NASA GCKP . computer 
program. This program calculates the time history of the 
initial chemical species concentrations, given the initial 
conditions of pressure, temperature and composition by 
integrating rate equations . Other constraints are also 
possible, such as keeping pressure and temperature unaltered. 
The input to this program requires , in addition to the 
initial conditions, a chemical reaction kinetic scheme to 
describe the combustion process. 

Depending on whether there is mixing of burnt gas and 
fresh air, to simulate addition of air, or not, i.e. 
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depending on whether the kinetic calculations are to be 
performed for the secondary and dilution zones, or the 
primary zone, respectively, the initial species concentra- 
tions are calculated in a slightly different manner. In 
the primary zone, the fuel and oxidizer are assumed to be 
perfectly mixed instantaneously and the initial mass 
fractions are calculated as follows . 

Dry Air ~ + ©2 + Ar + CO 2 

= 0.7555 + 0.2313 + 0.0127 + 0.0005 

= a. '0000 

Water Vapor = H„0 = Dry Air xW xRH = IxW xREJ 

^ s . s 

= W X RH . 

• s 

Fuel = CH^ = Dry Air x f/a = lx f/a = f/a . 

Total mass of initial species - (1 + W x RH + f/a) 

® ( 8 ) 

The mass fractions of each species ; 

1^2 = 0.7555 /(I + Wg X RH + f/a) , 

O 2 = 0.2313 /(I + Wg X RH + f/a) , 

Ar = 0.0127 /(I + W x RH + f/a) , 

s 

CO 2 = 0.0005 /(I + Wg X RH + f/a) , 

H^O = W X RH/(1 + W X RH + f/a) , and 

CH^ = (f/a) /(I + Wg X RH + f/a) . 
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A FORTRAN program given in Appendix 1 performs the above 
calculations and produces results on punched cards, ready 
for use in the chemical Jcinetic calculations computer 
program. 

When there is mixing of burnt gas and fresh air, as 
is the case during entry to the secondary and dilution 
zones, where the burnt gases from the previous zone are 
mixed with appropriate amounts of .air, to lower the fuel/air 
ratio, the initial species concentrations (mass fractions) 
for use in the kinetic calculations program, are calculated 
with the aid of the FORTRAN program in Appendix 2 as follows. 

Consider the case where it is required to calculate 
the initial species concentration during entry to the 
secondaary zone. Burnt gases from the primary zone are 
instantaneously mixed perfectly with fresh secondary air, 
and then allowed to burn in the secondary zone. The 
fuel/oxidizer ratios in the primary and secondary are 

(f/o) and (f/o) respectively. This situation can be 

P s 

considered as starting out with a certain amount of fuel, 
f, and the appropriate amount of oxidizer, namely f/(f/o)g, 
which will result in a fuel/oxidizer ratio of (f/o)^ in 
the secondary, feeding part of this oxidizer through the 
primary, namely f/Cf/o)^, and feeding in the rest of 
the oxidizer, namely f.{l/(f/o)g - l/Cf/o)^} as secondary 
air, which is mixed with the effluents from the primary 
zone combustion, to simulate entry into the secondary zone. 
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The mass fractions of the -burnt gas mixture are 
obtained from the results of plug flow calculations applied 
in the primary, zone. The mass fractions of the secondary 
air are the same as those for fresh air, depending on 

t 

ambient conditions. The mass of the burnt gas consists of 
the mass of primary air and injected fuel/- and the mass of 
the secondary air is the rest of the oxidizer mass as 
mentioned above. The masses of individual species in 
the burnt gas and secondary air are then given by multi- 
plying the total mass of either the burnt gas or secondary 
air, with their respective mass fractions. Instantaneous 
mixing ,of these two result in .a. different composition, 
calculated as given below. . 

Burnt gas individual species mass : 

rtu = H, X (f +'f/{f/oj ) (10) 

"i "“i . K 

Secondary oxidizer individual species mass : 

nr = H X (f/(f/o) - f/(f/o) ) (ID 

i i ^ 

After instantaneous mixing , individual species mass :■ 

nii = mj^_ + mo_ (12) 

After mixing total mass : 

M = I = I ' (13) 

i i i i 
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The mass fractions at secondary zone entry: 






I "Hd, 


+ m. 


+ m. 


H [1 + l/(f/o) ] + H^Jl/(f/o)g “ l/(f/o)p] 

^ i . 

g S- - ' ■' ' " 

^ ^ H. [1 + l/(f/o) 3 + H^_ll/{f/o)g - l/{f/o)p3 

( 14 ) 

These are the mass fractions at the entry to the 
secondary zone. The mass fractions at the entry to the 
dilution zone are calculated in the same manner ^ by con*- 
sidering the primary and secondary zones to be one unit, 
replacing the primary zone in the calculations-? and con- 
sidering the dilution z:one as the secondary zone. 

To describe the combustion process ? a kinetic chemical 
reaction scheme is also reguired. A list of reactions 
employed for this purpose is given in Table II. These 
reactions describe the guick breakdown of the methane fuel 
into radicals and carbon monoxide, and their subseguent 
oxidation to final products of water and carbon dioxide? 
and also other radical recombinations. In addition? the 
kinetics of the formation of oxides of nitrogen are also 
included. Table II also gives the forward rate constant 
associated with each reaction. Table Ila gives the 
different reaction schemes actually employed in the 
analytical program. 
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A 30 reaction scheme (number 1) was first employed 
in a pilot analytical effort to ascertain whether 
humidity had any effect on carbon monoxide production. 

This scheme is given and discussed in detail in a review 
paper by Osgerby^^^^ and includes reactions to describe 
burn up of methane, production and consumption of carbon 
monoxide , and production of both nitric oxide and 
nitrogen dioxide. 

A 20 reaction scheme (number 2 ) , which is essentially 

the same as the one used by Ay and Sichel^^*^^ , is also 

given in the same tables. The differences between 

schemes 1 and 2 are that the latter does not include the 

kinetics of nitrogen dioxide production, some auxilliary 

reactions for nitric oxide production and a carbon 

monoxide production reaction. Table Ila gives the details. 

In carbon monoxide kinetics , the reaction CO + OH = CO 2 + H 

has been shown to dominate the oxidation -process by 
(25) 

Westenberg , and nitric oxide kinetics have been 

adequately described by the mechanism proposed by Zeldovich ^ 

The reaction scheme 3 consists of the same 20 

reactions , but with a modified forward rate constant for the 

CO + OH = CO 2 + H reaction. This rate constant is similar 

(27) 

to the one developed by Kollrack ^ ' and was chosen for 

its large temperature dependence . 
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A. Analytical Models t 

A typical gas turbine coinbustor is shown in Fig. 1, 
wherein the combustor is subdivided into the three familiar 
zones, the primary, secondary and dilution. Combustion 
occurs in each of these three zones accompanied by complex 
flow processes. In the analytical model, these processes 
are ignored and simple plug- flow reactors are employed 
to simulate combustion in the three zones, one for 

each zone as shotm in Pig. 8. 

In the primary zone , fuel and oxidizer , taken in 

appropriate amounts to simulate a. particular fuel/air 

ratio in the primary zone, are instantaneously mixed and 

allowed to undergo chemical reaction at a pressure of 

4 atmospheres (compressor pressure .ratio is four) and at 

a temperature corresponding to the adiabatic ‘flame 

temperature which reflects both ambient conditions and 

fuel/air ratio. • The reaction process is allowed to proceed 

for 5 msec., which is taken as representative of a primary 

zone residence time. The initial species composition is 

calculated using the FORTRAH program given in Appendix 1, 

and the chemical composition of the primary zone effluents, 

using the NASA GCKP computer program, developed by 

(23) 

Bittker and Scullin 

These effluents are then perfectly mixed, instan- 
taneously, with precalculated amounts of fresh air in 
the manner explained in the previous section, using the 
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PORTRMf program given in Appendix 2 'to simulate -entry 
into the secondary zone. The fuel/air ratio in the 
secondary zone is lower than that in the ' primary zone and 
hence the temperature is also lower. Using the same 
NASA GCKP program, the chemical composition of the mixture 
after a secondary zone residence time of 5 -msec, is 
calculated. The conditions in the secondary zone ‘are 
once again modeled as' a plug-flow reactor. 'The secondary 
zone effluents are once again mixed with an additional 
amomt of fresh air to simulate entsry into the dilution 
zone. in the same manner as above, and the mixture' allowed 
to undergo chemical reaction in the plug- flow reactor 
mode-ling, the dilution zone, and the chemical composition 
calculated after a dilution zone -residence time of 5 ‘msec. 
This composition is taken to be that -at the combustor exit. 
Carbon monoxide mole fractions, reflecting various 
ambient conditions' -and fuel/air ratios, in particular, 
are -noted. 

■■ This type of combustor modelling applies for schemes 
which assume' the primary -zone to be fuel-rich. However, 
as stoichiometric primary zone fuel/air ratios are 
approached, a slightly modified combustor model is necessary. 
The details of the modif ications are given in what 
follows . 
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It is suggested that at near and below stoichiometric 

primary zone fuel/air ratios, homogeneous single phase 

combustion may not occur, but that individual fuel droplet 

burning may dominate the combustion. Actual burning 

occurs around individual droplets under stoichiometric 

conditions and the products of this combustion are 

immediately quenched by the surrounding air to lower 

fuel/air ratios, in the primary zone itself. This has the 

effect of "freezing" the highest carbon monoxide levels 

which occur during droplet combustion. Indeed, experi- 

f 28 ) 

mental evidence given by Sullivan , indicates that this 
type of primary zone quenching does occur. Sullivan found 
that the exact placement of air holes in the secondary zone had 
very little effect on the carbon monoxide emission levels, 
suggesting that these levels are quenched by surrounding 
primary air itself and any additional air would make very- 
little difference except in lowering the mass fractions 
because an additional mass of air is introduced. 

To model this type of behavior, the primary zone is 
considered as comprising two combustion phases. The first 
phase is the individual droplet burning and the second 
phase is the immediare dilution by surrounding primary air. 

As shown in Fig. 9, each such zone is modeled as a plug 
flow reactor. The first reactor has a temperature and 
initial species concentration reflecting ambient conditions 
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and for a stoichiometric or slightly rich fuel/air ratio. 
Chemical reaction occurs in this zone for a very short 
period of time, 0.1 msec., which is of the same order 
as the time taken by a flame front to travel through the 
near~stoichiometric combustion zone surrounding 
individual droplets. The products of this combustion 
zone are now mixed with fresh air to bring the fuel/air 
ratio, and hence the temperature down, simulating the 
process of immediate dilution of droplet buiming products, 
and the mixtxxre allowed to undergo chemical reaction in 
the second primary plug flow reactor for' 5 msec . , 
representative of a primary zone residence time. 

Subsequently, these products are allowed to pass 
through the secondary and dilution zones again modeled 
as plug flow reactors , in the same manner as for the 
previous model. .The chemical composition, carbon monoxide 
mole fraction in particular, is noted at the end of the 
dilution zone. 

In the analytical model presented here, various fuel/air 
ratio combinations were considered, and the various para- 
meters relevant to one of these values, namely an ^ of 
0.011, are given in Table III. Experimentally, the combustor 
was operated under conditions corresponding to the primary, 
secondary, and dilution zone fuel/air ratios 'given in the 
table. For an overall or dilution zone fuel/air ratio of 
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0.0-11/ the secondary and primary zone fuel/air ratios 
lie between 0.034 and 0.019/ and' 0.071 and 0.048, 
respectively. A comparison of the analytical and typical 
values of fuel/air ratio- reveals a sufficient. range in 
the variation of the typical values to include analytical 
values. However, differences in the residence times in 
each zone are -apparent. In the primary zone, the actual 
residence is between 1.90 and 1.60 msec, whereas the 
analytical model takes this as 5.0 msec. This may be 
.justifiable because chemical equilibrium is quickly 
attained in the primary zone , due to large reaction rates , 
and there is practically no difference between the carbon 
monoxide mole' fraction after 1-0 msec, or after 5.0 msec, 
residence time in the primary .zone, as per calculations. 
This situation is illustrated in Pig. 10, wtiichTTs a plot 
of the carbon monoxide mole fraction pbtained from ■ 
eguilibrixmi ' and, kinetic calculation for various fuel/air 
ratios. In the entire range of values investigated, there 
is 'very little difference between the two unoJLe fractions .- 
However, deviations become apparent as .either • very .lean or 
very rich mixtures are considered. 

The secondary zone residence time can be seen from 
Table III, to be- of the same order in both .the typical 
combustor and the analytical model. It will be shown later 
in the section on analytical results, that the secondary 
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zone parameters play a dominant role in determining the 
carbon monoxide' emission levels and hence the closeness 
in both fuel/ air ratio and residence time in this zone 
is, indeed, necessary. 

In the dilution zone, the residence time taken in 
the analytical model is about twice the typical value. 

The temperatures in the dilution zone are relatively low , 
with correspondingly lower reaction rates , leading to 
"frozen” chemistry. There is thus little change in the 
carbon monoxide mole fraction between 2,. 5 msec, and 5.0 
msec, dilution zone residence time and hence the choice of 
the larger residence time may be justifiable. 

Further, in the analytical model, methane was used 
as the fuel instead of JP4 , which was used in the .experi- 
mental program. It has already been pointed out that the 
primary zone carbon monoxide levels are close ±o equilibrium 
levels, as apparent from Fig. 10. It may be seen from 
Table IV that the equilibrium carbon monoxide mole fractions, 
for the same conditions, employing methane and JP4 as fuel, 
are of the same order, but JP4 results in a slightly higher 
(about one-and-a-half times) mole fraction. However, the 
adiabatic flame temperatures for methane are slightly 
lower than those for JP4, for the same conditions, 
resulting in slower carbon monoxide burnup when methane 
kinetics are employed. It is possible that these two 
effects may compensate for each other, at least to a limited 
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extent, resulting in more realistic carbon monoxide 
levels. It is recognized that using methane instead of 
JP4'is, indeed, an approximation, but because the kinetics 
of JP4 -oxidation are only poorly understood, use of 
methane as the fuel may also be justifiable. 

As rega-rds' the chemical kinetics employed in the 

analytical effort, three schemes have already, been 

presented. The scheme nvimber 1 has been found to be 

imsatis factory in predicting carbon monoxide levels, in 

that the calculated mole fractions were two orders of 

magnitude smaller than measured values in some cases . 

Further ,• scheme number 1 also contained reactions describing 

the production of nitrogen dioxide, which ;was found to 

be several orders of magnitude smaller than nitric oxide , 

as a result of calculations . Reactions involving nitrogen 

dioxide were therefore considered 'unimportant and reaction 

scheme number 2 , consisting of 20 reactions given by Ay and 
(2A) 

Sichel , which did not contain .such reactions was then 
used. However, its use was limited to the primary zone 
alone, since it too underpredicted carbon monoxide mole 
fraction at the end of the dilution zone ,by two orders of 
ma^itude. Further, the dilution zone carbon monoxide 
levels did not exhibit increasing trends with humidity, 
when this scheme was used. This led to a search for a 
rate constant for the carbon monoxide oxidation reaction 
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CO + OH = CO^ + H, which is highly temperature dependent, 
because ambient humidity causes a decrease in the flame 
temperature and it is required to reduce burn-up of 
carbon monoxide sufficiently, as humidity is increased, 
in order to predict measured trends. The search resulted, 
in the use of a rate constant developed by Kollrack . 

This rate constant had the seventh power of the absolute 
temperature in the pre-exponent factor, indicating high 
temperature .sensitivity. However, when thiS' rate constant 
was used in the analysis, the carbon monoxide emission 
•levels calculated were once again two orders of magnitude 
smaller than the measured levels. It was then suggested 
that this rate constant be reduced by a factor of ten,- and 
when this value was used,’ reasonable results were obtained. 
Figure 11 illustrates this situation adequately"^ For a 
• variety of conditions , the carbon monoxide emission level 
is -plotted against the combustor exit temperature (dilution 
zone flame temperature) . Comparison of measured carbon 
monoxide ppm levels and calculated values reveals not only 
good agreement in the order of magnitude, but also in the 
trends as regards the effects of ambient himidity and 
temperature . 

Ftirther, Table V gives the values of the largest and 
the smallest reported , rate constants for the carbon 
monoxide oxidation reaction for various absolute temperatures . 
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These rate constants are taken from the work of Sorenson 
et al. ' . In the secondary zone, the temperature is of 
the order of 1500 °K, at which the used rate constant 
"USED K" is slightly smaller than the smallest reported 
value, "SMALL K" . However in the dilution zone, the 
temperature is very low, of the order of 800 ®K, at which 
the used rate constant is slightly larger than the 
smallest value but smaller than the largest value, "LARGE K" , 
The used rate constant is, indeed, relatively small, 
but because of good results as in Fig. 11, its use may 
also be justifiable. 

B. Analytical Results : 

The experimental results, already being ascertained, 
in that the principal effects of ambient temperature and 
humidity on carbon monoxide emission levels at the com- 
•b us tor exit are as follows : 

(1) ambient temperature increases result in decreases 
in carbon monoxide level, for zero hiimidity and 
for a particular overall fuel/air ratio and 
compressor discharge condition, and 

(2) ambient humidity increases cause increases in 
carbon monoxide emission levels, for a particular 
ambient temperature, overall fuel/air ratio, and 
compressor discharge conditions . 
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It has been possible to predict the above trends 
employing the analytical model described in the previous 
sscti'On , to a reasonable extent. The calculated results 
are presented here for a compressor pressure ratio of 
four. 

Results of the adiabatic flame temperature calculations 
using the NASA CEC71 ' computer program for various 
equivalence ratios from 0.2 to 1.3, for ambient temperatures 
of 244, '289 and 322°K and for relative humidities of 0, ,50 
and 100 percent are presented in Fig. 12.. The well-known 
effects of ambient temperature and humidity on the 
adiabatic flame temperature are apparent. Ambient tern- ' 
perature increases produce increases in the flame tempera- 
ture, for a particular equivalence .ratio and zero humidity. 
Ambient humidity increases cause decreases in the flame 
temperature at any particular equivalence ratio and ambient 
temperature. For an ambient temperature of -244®K, Pig. -12 
gives only one curve; this is because at this ambient 
temperature, a very small amount of water vapor corresponds 
to saturation and there is therefore, practically no effect 
on the flame temperature. 

The effect of humidity on fuel/air ratios corresponding 
to particular equivalence ratios for various' ambient 
humidities is given in Fig . 13 . Near the stoichiometric 
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equivalence ratio this effect is not seen, but as .one 
deviates from stoichiometric values, the effect is quite 
pronounced in that below stoichiometric values , the 
fuel/air ratio for humid mixtures is less than that for 
dry mixtures and above stoichiometric values , the 
reverse is true. Figure 13 contains the fuel/air values 
plotted against equivalence ratio for ambient temperatures 
of 244, 289 and 322 ‘’K and for relative humidities of 
0, 50 and 100 percent for the same reason as mentioned 
above , there is only one curve for an ambient temperature 
of 244°K. The curves in Fig. 13 are useful when it is 
desired to determine flame temperatures for a particular 
fuel/air ratio, instead of for a particular equivalence 
ratio. Since flame temperatures in Fig. 12 are plotted 
against equivalence ratio. Fig. 13 can be used to determine 
the equivalence ratio values corresponding to a particular 
fuel/air ratio and various desired ambient conditions of 
temperature and humidity. Using these equivalence ratio 
values and corresponding curves in Fig. 12, the flame 
temperatures may be determined for the various conditions 
considered. The combustor is modeled as three plug-flow 
zones, primary, secondary and dilution, and results from 
calculations for each zone and relevant discussion follow. 

The primary zone results are given in Fig. 14, which 
is a plot of the carbon monoxide correction f actor. CF__, 
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against the ambient temperature', for relative humidities 
of 0, 50 and 100 percent, and for primary zone fuel/air 
ratios of 0.08, 0.07, and 0.06, the last of which is 
nearly stoichiometric.' The' correction factor has been 
defined as : 

QY Carbon monoxide mole fraction at standard conditions 

CO .Carbon monoxide mole fraction at nonstandard conditions 

with standard conditions taken as 289 °K temperature and 
zero humidity. The mole fractions are taken after allowing 
fuel and oxidizer to undergo chemical reaction in the 
primary zone plug flow reactor for a residence time of . 

5 msec. 

The three sets of curves in Fig. 14, one for each 
primary zone fuel/air xatio, are similar ,in nature, except 
for differences in the slopes. For zero h'umidity, in all 
three sets of curves, the slope of the CF^^ curve is 
negative, indicating that the carbon monoxide levels 
increase with .increasing ambient temperature at zero 
humidity. At the two other relative humidities, 50 and 
100%, ‘the CFj^q curves have positive .slopes, indicating that 
at these humidities , an increase in ambient temperature 
leads to a decrease in the carbon monoxide level. At 
any particular temperature, increasing humidity leads 'to 
an increase in CF^q, indicating that the carbon monoxide 
level decreases as himidity is increased at a particular 
ambient temperature . 
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Further, products of the primary zone plug-flow re- 
actor are almost in equilibriim, as has been mentioned 
earlier, with reference to Fig. 10, and thus both 
equilibrium and kinetic results indicate carbon monoxide 
trends which are to increase with ambient temperature 
-and to decrease with ambient humidity. These trends are 
precisely the inverse of those measured at the combustor 
exit, but have been experimentally verified by Muller- 
Dethlefs -and Schlader^^^^ , who made measurements in rich 
burner flames. which correspond to primary zones. These 
results may be explained by considering the effect of 
flame temperature on dissociation. It is therefore, 
logical to assume, in view of these results, that the 
subsequent secondary and dilution zones are instrumental 
in reversing these trends. Further, it is evident from 
a perusal of Fig. 14, that there are larger deviations 
near stoichiometric fuel/air ratios, making it increasingly 
difficult for the secondary and dilution zones to reverse 
the trends . This explains why the rich primary zone 
model (Pig. 8) shows signs of failure for near-stoichiometric 
primary zones . 

The temperatures in the dilution zone are lower, hence 
chemical composition is somewhat "frozen"; so this zone may 
not have a great effect in, reversing the trends. The only 
zone left is the secondary zone, where the temperatures 
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are not as low. Indeed, it turns out that the secondary 
zone can have a dramatic effect on the carbon monoxide 
levels. This is illustrated in Fig. 15, which is a' plot 
of the carbon monoxide emission level in mole fractions 
at the combustor exit, against the exit temperature 
calculated for various secondary ftiel/air ratio com- ■ 
binations and for various ambient conditions as shown in 
-the figure, employing the analytical model for rich 
primary zones (Fig. 8) . The primary zone and dilution 
zone fuel/air ratios are fixed, 0.07 and 0.015 respectively, 
but the intermediate secondary zone fuel/aif ratio is 
varied -as follows: 0.0 35 , 0.03, 0.0275 and 0.025., 

However, when secondary fuel/air ratios beyond the range 
between 0.0 35 and 0.0 25,, namely 0.04 and 0.0 2 for instance, 
are considered, the results fail -to predict the trends 
mentioned above. The results are plotted in groups of 
symbols as explained in the 'figure , each group corresponding 
to one analytical scheme comprising the three combustion 
zones with residence times of 5 msec in each zone. Changing 
the secondary zone fuel/air ratio from 0.035 to 0.025 re- 
sults in an order of magnitude change in the carbon 
monoxide mole fraction. It is thus confirmed that the 
secondary zone is the controlling factor as regards the 
carbon monoxide emission level. 



While the effect of the secondary zone is evident, 
it is also noticed that the carbon monoxide trends with 
respect to ambient temperature and humidity have also 
been reversed. The results on Fig. 15 indicate that for 
each of these fuel/air ratio combinations, the analytical 
.model is capable of predicting the observed carbon monoxide 
emission trends, in that for a particular analytical 
scheme, the carbon monoxide decreases with increasing 
ambient temperature (compare identically shaded circular 
symbols) and increases with increasing hxmidity at a 
particular ambient temperature (compare identically shaded 
and, flagged symbols) . Further, it may be noted that the 
largest levels of carbon monoxide are produced for .smaller 
values of secondary fuel/air ratio, which occurs because 
the high CO levels in the primary zone are "frozen" by 
these 'Secondary zones , where the lower temperature 
"quenches" the CO oxidation reaction. 

Employing the same analytical model for rich primary 
zones (Fig. 8) , additional calculations over a wide range 
of conditions yielded the same kind of results. Figures 
16, 17 and 18 are plots of carbon monoxide mole fraction 
against the combustor exit temperature, which is the dilution 
zo.ne adiabatic flame temperature. Three groups of results 
are presented in each figure, one for each fuel/air ratio 
combination. ’ All three figures start out with rich primary 
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zones, considerably different from stoichiometric values, 
and subsequently undergo secondary and dilution zone 
combustion as follows : 


^/a)p 

f/a)^ 

f/a)g 

0 .070 

0.030 

0.020 

0..070 

0.0275 

0.015 

0.070 

0 .015 

0.010 

0.080 

0.035 

0.020 

0 .080 

0.030 

0.015 

0.080 • 

0.025 

0.010 

0.090 

0.035 

0.020 

0.090 

0.0 30 

0.015 

0.090 

0.025 

0.010 


The residence times in each zone was taken as 5 msec, and 
calculations were performed for ambient conditions of 
244, 289 and 322°K temperature and 0, 50 and 100% relative 
humidity . In each individual case , other secondary 
fuel/air ratio combinations were also tried, but the 
analytical model gave reasonable results only for a 
limited range of secondary fuel/air ratio . . -Only such 
values of secondary fuel/air ratio and' accompanying primary 
and dilution fuel/ air ratio combinations are presented 
here. Further, these values correspond closely with 
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with typical combustor values , where the combustor fuel/ . 
air ratio ranged from 0.015 to 0.007. 

On each figure, the effects of ambient conditions 
are apparent in that the trends obtained from calculations 
agree' with experimentally observed trends as regards 
carbon monoxide emission levels. However, there are 
differences in the degree of these effects, depending on 
the .particular choice of fuel/air ratios in the three 
plug flow zones . -Generally , richer primary zones may be 
associated with richer secondary and dilution zones; 
doing so, .and comparing results from the three figures 
16, ,17 and 18 for the following coiribinations : 


f/a)p 

f/a)^ 

f/a)a 

0.070 

0.025 

0.010 

0.0 80 

0..030 

0 .0 15 

0.090 

0.035 

0.020 


It may be seen that rich primary .zones coupled with 
richer secondary and dilution zones result in lower carbon 
monoxide emissions in general. Also, the slopes 

(p= — ) and .(-Tm — ) are seen to depend on the 

'^4 RH^O.O -4 HH=1.0 

fuel/air ratio combination. For the same primary zone 

fuel/air ratio', and for higher values of secondary and 

dilution zone fuel/air ratio, both these slopes can be 

seen to decrease. 
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It was previously mentioned, in the section on the 
analytical models , that the model employed was only 
successful as applied to sequences starting out with rich 
primary zones . Lean primary zones may be handled in a 
similar manner, but with the modification that individual 
droplets of fuel first burn in the primary and are 
immediately diluted down to lean fuel/air ratios. This 
scheme, for lean primary zones is given in Fig. 9 and has 
been discussed already. 

Employing this scheme then, results for two sequences 

X 

starting with a primary zone with a fuel/air ratio of 
0.06, which is close to the- stoichiometric value, and 
subsequent dilution to a fuel/air ratio -of 0-03 in the 
primary zone itself and then completing the combustion 
process through secondary and dilution zones with fuel/air 
ratios of 0.02, 0.015 and 0.010, 0.007 respectively, are 
shown in Fig. 19., That this model is capable of predicting 
the same CO trends with regard to ambient temperature and 

hxmaidity changes , is evident in this figure . 

^ \ 

'Combined results for both the analytical models are 
given in Fig. 20, which has the carbon monoxide mole 
fraction plotted against dilution zone temperature, for 
the following sequences: 
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f/a) 

Pi 

f/a) 

P2 - 

f/a)^ - 

f/a) 3 

0.090 

— 

0.035 

0 .015 

0.080 

— 

0.030- 

0.010 

0.060 

0.030 

0.015. 

0.007 


The points in the figure are in, groups^ one for each 
dilution sequence and the effects of ambient temperature 
and hiJmidity can be seen to be similar for each sequence 
and agree with measurements.,- Care has been taken to 
couple richer primary zones with richer subsequent zones 
and vice versa. Once again, richer zones produce less 
carbon monoxide levels. These combinations fit in .well 
with -experimentally determined values of fuel/air .ratio. 

The analytical models have .been applied over a wide 
range of operating conditions and have been reasonably 
successful in predicting carbon monoxide emission levels 
both as far as order of magnitude and trends with .ambient, 
conditions are concerned. In a manner s'imilar to the 
experimental results , these analytical results may be 
summarized as follows: 

As regards the carbon monoxide (mole fraction) 
emission level: 

(i) For typical combustor conditions represented in the 
analytical model, the calculated- carbon monoxide mole 
. fraction (converted to ppm) and the measured values are 


44 



roughly the same because the rate constant for the carbon 
monoxide oxidation reaction was modified, to produce, 
these results . 

(ii) For a fixed set of ambient conditions, increasing 
the secondary zone fuel/air ratio (and hence the dilution 
zone fuel/ air ratio) results in a decrease in CO level. 

(iii) For a fixed fuel/air ratio sequence, and zero 
humidity, increasing the ambient temperature causes the 
CO mole fraction at combustor exit to decrease. 

(iv) For a fixed fuel/air ratio sequence and ambient 
temperature, an increase in the humidity leads to an 
increase in the CO emission level. 


(V) 


ax 


The slope ( 


CO 


3T. 


increases with decreasing 


RH=0.. 0 

dilution zone fuel/air ratio for any (f/a) and (f/a) . 

p s 

ax. 


(vi) The slope ( 


CO 


r>m / generally increases with 

^^4 RH=1'.0 


decreasing dilution zone fuel/air ratio, for any .(f/a) 


and (f/a) 
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ANALYTICAL AND EXPERIMENTAL COMPARISONS 

From the results of the experimental and analytical 
program, it may be’ seen that, indeed, plug-flow reactor 
modelling of the gas turbine combustor can be reasonably 
successful in predicting both carbon monoxide emission 
levels - with the aid of a reduced rate constant for the 
carbon monoxide oxidation' reaction - and the effects of 
ambient temperature and humidity on carbon monoxide 
emission, with the aid of the large temperature dependence 
in the same rate constant. 

Figure 11 presents the experimentally measured and 
calculated carbon monoxide emission levels (ppm) for similar 
combustor discharge conditions . It can be seen that 
the amounts of carbon monoxide calculated and measured are 
roughly the same , over a fairly wide range of conditions . 

Figures 21, 22 and 23 present measured and calculated 
carbon monoxide emission levels in terms of the correction 
factor against ambient temperature for relative humidities 
of 0, 50 and 100 percent, and for various overall fuel/air 
ratios, which correspond to the dilution zone fuel/air 
ratio in the calculations. It must be kept in mind, 
howe'ver, that in the analytical model, each dilution zone 
fuel/ air ratio is associated with secondary and primary 
zone fuel/air ratios also. These values are given in 
•the figures. 
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In all these figures/ the similarity in the C_ 

. ^CO 

curves is striking. At zero humidity, these curves have 
a positive slope, in all three figures/ indicating -chat 
the carbon monoxide- levels decrease with increasing 
ambient temperature. The slopes of these curves for both 
analytical and experimental results are about the same. 

However, at relative humidities other than zero, 

(for 50 and 100 percent) , the curves obtained from the 
.analytical program do not indicate as large an increase in 
carbon monoxide production as produced during the 
experimental .measurements . The analytical model somewhat 
.underpredicts the effect of humidity on carbon monoxide, 
but is, at least, successful .in predicting .similar trends. 

Results from other experimental programs , such as 
those presented by .Blazowski et al . include carbon 
monoxide correction factors for ambient temperature for 
two different fuels JP^ and JPg and for two design 'varia- 
tions of a T-56 combustor, employing a stoichiometric -and 
a lean primary zone. These results are plotted in Figure 
24 as the correction factor against ambient temperature. 
"Least squares" curve fits to these data, employing 
linear interpolation resulted in equations for these 
correction factors, given in Table VI. The curve fitting 
was done using the BMD programs package . The equations 
give the correction factor in terms of. the pressure ratio 
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and ambient, temperature. The effect of humidity was not 
investigated here, but some confidence can be gained by 
comparing the slopes of the zero humidity correction factor 
curves in Figures 21, 22 and 23 with the curves in Figure 
24. The slopes and the shapes of the curves are surprisingly 
similar. 
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CONCLUSIONS 


From the foregoing discussion of the experimental 
and analytical results , the following conclusions may 
be made: 

1. The. analytical models presented can, with limited 
success, predict carbon monoxide emission levels from 
idling gas turbine combustors. Predicted and measured 
levels are roughly the same . This agreement was forced 

by modifying the rate constant for the carbon monoxide 
oxidation reaction. 

2 . As regards prediction of the -effects of ambient 
temperature and .humidity on carbon monoxide emissions, 

the analytical model is able to predict the measured trends , 
but is unable to exactly match the slopes of these trends, 
in the sense, the model underpredicts the -'effect of 
humidity on carbon monoxide emissions . The effects of 
temperature, however, are in reasonable agreement. 

3. For zero humidity, the correction factors for 
both the JT8D~17 and the T-56 combustor are similar; so it 
may be possible to have a general equation for correction 
factors for ambient temperature for both engines . It is ■ ’ 
one of the goals of this- research effort to examine such a 

/possibility, to help ease the regulatory task of developing 
corrections for nonstandard inlet conditions . 

I 
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4. With simple plug-flow combustor modelling 
these results have been possible - one of the goals of 
this research effort was simplicity, which has been 
retained in the modelling. 

5. The results indicate that carbon monoxide pro- 
duction in quantity is due to quenching of the primary 
zone products by lean secondary zones , resulting in 
"super-equilibrium" carbon monoxide emission levels. 
Gentle dilution of the primary zone will help bum-up 

of the carbon monoxide to a larger extent and thus result 
in lowered emissions . This is in agreement with results 

( 31 \ 

of Palmer and Seery 
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TABLE I 


idle JT8D - 17 COMBOSTOR CONDITIONS 

Nominal Operation: 

Total Inlet Pressure 
Total Inlet Temperature 
Air Flow 
Fuel Flow 
Fuel/Air Ratio 

Test Conditions; 

Compressor Efficiency - 0.8 

Compressor Pressure Ratio 2, 3, 4, 5 

Compressor Inlet Pressure 1 atm. 
Compressor Inlet Temperature 244, 289, 322 K 
Compressor Inlet Relative 


• 

Humidity 

0, 50, 

100 % 

Fuel/Air 

Ratio 

.015, 

. OIL, 

Constant 

Compressor Discharge 




Mach Niunber 

0.42, 

or, 


Constant Reference Velocity 15,2 m/sec 


2,47 atm. 

393 K 
1.37, kg/sec 
.0161 kg/ sec 
.0117 



TA3L3 HI 


Kinetic Scheme for Methane/Air Combustion With 




Forward 

Rate Constants : 

k- = A e 

-AE/RT 











units : cm. 

y mo 1 e / 

sec. 

No. 



R€ 

;acuion ■ 



A 

a 

■ AE 

1 

M 


CH^ 

= 

CH 3 


H *f 

14 ’ • 

0.20E18 

0.0 

88332.5 

2 

CH, 


'OE 

— 

CH 3 


H 2 O 


0.28E14 

0.0 

4962.5 

3 



0 

- 

CE 3 

-r 

OH 


0 .20E1.4 

0.0 

9210.4 

4 


-f 

K 

= 

CH 3 

-r 

^^2 


0.69E14 

0.0 

11810.8 

5 

CH 3 

-f 

^2 

c; 

HCO 

-f 

H 2 O 


0.20E11 

0.0 

0.0 

6 

CE 3 

H- 

-0 

= 

HCO 


^^2 


0 . 10E15 

0.0 

0.0 

7 

HCO 


OE 

= 

CO' 

- 7 - 

H 2 O 


0 .10E15 

0.0 

0.0 

8 


-f 

HCO 

= 

CO 

-r 

H + 

M 

' 0.20E13 

0.5 

28584.0 

9 

CO 

JL 

OH 

- 

C 02 


H 


‘0 . 56E12 

0.0 

600.0 










0 . 85E-14 

7.0 

-13895.0 

10 

E 


°2 

= 

0 


OE 


0.22E15. 

0 ..0 

16554.9 

11 

0 

-r 

^2 


H 

- 

OE 


0 . 1-7E14 

0.0 

9428.8 

12 

0 


E 2 O 

- 

OH 


OH . 


0.58E14- 

■0.0 

18004.0 

13 

B 

•f 

■^20 

= 

^-2 

*f 

OH 


0-.-84E14 

0.0 ' 

20048.5 

14 

E 

“f 

OH 

-f 

JE4 


E 2 O -4 

N 

0.40E20 

-1.0 

0.0 

15 

H 

-f 

E 

u. 

M 


E 2 4 

N 

0.. 15E19 

-1.0 

0.0 

16 

0 

■f 

0 

.j. 

M 

:= 

O 2 4 

M 

0.40E18 

-1.0 

0..0 

17 

0 


H 

+ 

M 

- 

OH 4 

M 

0.53E16 

0.0 

'5516. 3 

18 . 

N 

•f 

^2 

= 

NO 

*f 

0 


0.64E10 

1.0 

6232.9 

19 

. 0 

+ 

^2 

- 

NO 


N 


0 .14E15 

0.0 

'75231.5 

20 

OH 

JL 

N 


NO 

4- 

H 


0.40E14 

0.0 

0.0 

21 

0 


CO 2 

= 

CO 

•r 

°2 


0.19E14 

0.0 

54150 .'0 



TABLE ±1 (contirmed) 


No. 



He 

^action 




22 

^^2 

*f 

^2 


N 

-f 

0 



23 

^'2 


°2 

= 

NO 


NO 



24 

NO 

+ 

NO 

= 

N 

-f 

NO 2 



25 

M 


NO 

- 

N 

T* 

0 

•r 

M 

26 

0 

-r 

NO 

-~ 

Ki 

= 

NO 2 


M 

27 

M 

*f 

NO 2 

= 

N 

-r 

°2 

-r 

N 

28 

0 


NO 2 

- 

NO 





29 

H 

-i- 

0 

to 

- 

NO 


OE 



30 

N 


C 02 

= 

CO 


NO 



31 

CO 

•f 

N02 

= 

NO 

•f 

CO 2 




A 

a 

AE 

0.27E15 

0 

1 

120428.0 

0.42E15 

0.0 

119100.0 

0.30S12 

0.0 

0 

0 

2.27E17 

-0 .5 

148846.0 

1.05E15 

0.0 

- -1870.0 

6 .00E14 

-1.5 

105200.0 

1.00E13 

0.0 

600.0 

7.25E14 

0.0 

19 30.0 

2.00S11 

-0.5 

7950.0 

,2.00E11 

0 

l 

4968.0 



TABLE II a 


REACTION SCHEMES 


Scheme: 

1 


2 


3 


R.eactioiiS : 


All reactions in .Table II 
except reaction 6 . 


The first 20 reactions in 
Table II, with the first 
given rate constant- for 
reaction 9 . 

The first 20 reactions in 
Table II, with -the second 
given rate constant for 
reaction 9.. 



TABLE III 


LOCAL PUEL/AIR RATIOS AND RESIDENCE TIMES 


TYPICAL COMBUSTOR: 



Primary 

Secondary 

Dilution 


High Low 

High 

Low 

High 

Low 

Fuel/Air Ratio 

.071 .048 

.034 

.019 

.012 

.011 

Residence Time 
{milliseconds ) 

1.90 1.60 

4.47 

3.54 

2.62 

■ 2.45 

ANALYTICAL MODEL; 






Primary 

Secondary 

Dilution 



High 

Low 

High 

Low - 

Pue 1/Air Ratio 

.07 

.035 

.025 

.020 

.010 

Residence Time 
(milliseconds) 

5.0 

5.0 

5.0 

5.0 

5.0 



TABLE IV 


EQUILIBRIUM CARBON MONOXIDE MOLE FRACTION AND ADIABATIC FLAME TEMPERATURE DEG 

PR = 4 ’METHANE FUEL 


Equivalence Ratio 

1.2 

1 . d ■ 

0.8 

1.2 

1.0 

0.8 

T K 
o 

RH % 

•' Equilibrium Carbon 

Monojcide 

Adiabatic 

Flame 

Temperature 

244 

0 

.04548 

.00795 

.00043 

2196 

2298 

2059 

244 

50 , 

.04548 

.00794 

.00042 

2196 

2298 

2059 

244 

ioo 

.04547 ' 

. 00793 

.00042 

2196 

2298 

2058 

289 

0 

.04584 

,; 00904 

.00060 

2238 

2331 

2103 

289 

50 

.04563 

. 00854 

.00050 

2221 

2319 

2083 

289 

100 

.04541 ■ 

. .00807 

■ . 00042 

2204 

2307 

2064 

322 

0 

.06169 

i 00991 

; 00076 

2269 

2,355 

2136 

322 

50 

.04448 

i 00662 

: 00022 

2151 

2269 

1998 

322 

100 

.04268 

.00440 . 

, .00006 

2045 

2190 

1872 



TABLE IV Conb'd 


EQUILIBRIUM CARBON MONOXIDE MOLE FRACTION AND ADIABATIC FLAME TEMPERATURE 

PR = 4 JP 4 FUEL 


Equivalence Ratio 

1-1 

1.0 

0.8 

1.2 

1.0 

o 

• 

CO 

T K 
o 

RH % 

Equilibrium Carbon 

Monoxide 

Adiabatic 

Flame 

Temperature 

244 

0 

.05948 

i 01284 

;00091 

2280 

2354 

2118 

244 

50 

.05948 

.01283 

. 00090 

2280 

2354 

2118 

244 

100 

.05948 

.01281 

.00090 

2279 

2354 

2117 

289 

0 

,05988 

. 01442 

.00125 

2322 

2386 

2162 

289 

50 

.05972 

;01367 

.00106 

2304 

2373 

2141 

289 

100 

.05955 

.01295 

.00089 

2287 

2361 

2121 

322 

0 

.06022 

-01567 

.00156 

2353 

2409 

2194 

322 

50 

;05886 

.01069 

. 00047 

2229 

2322 

2050 

322 

100 

.05722 

.00725 

'. 00013 

2117 

2241 

1918 



TABLE V 


RATE CONSTANTS FOR THE 

-,CO + QH =.C02 + H 

REACTION 

T,Deg. 

Large K 

Used K 

Small K 

Kelvin 




500.00 

.22833E10 

.79860E11 

. 118-9 2E10 

600.00 

.12863E11 

.27749E11 

.30957E10 

700.00 

.44221E11 

,1S419E11 

.61311E10 

800.00 

.11165E12 

.11249E11 

.10236E11 

900. 00 

.22944E12 

.97043E10 

.15249E11 

1000.0 

.40826E12 

•.93214E10 

.20977E11 

1100. .0 

.65419E12 

. 96130E10 

..2723 0E11 

1200.0 

.96903E12 

.10401E11 

.33844E11 

1300.0 . 

.13512E13 

.11629E11 

.40680E11 

1400.0 

.17967E13 

.13298E11 

.47629E11 

1500.0 

.23000E13 

.15444E11 

.54604E11 

1600.0 

.28548E13 

. 181-2 6E11 

.6154IE11 

1700.0 

.34546E13 

.21421E11 

.68389E11 



TABLE VI 


CORRECTION FACTORS, CARBON MONOXIDE, T 56 COMBUSTOR 
PR = Pressure Ratio T^ = Ambient Temperatur 

Primary Zone Fuel Correction Factor 

Std. JP 4 .970772 ~ .045481*PR 

+ .001057*PR*T 

o 


std. 


JP 8 1.031017 


.05713*PR 

.000843*PR*T 


Lean 


JP 4 


1.030849 - 
+ 


. 0S9991*PR 

. 000961*PR*T^ 
o 


Lean 


JP 8 


1.027965 - 
+ 


. 050659*PR 
,000737*PR*T 

o 



appendix 1 


program FDR CQHPUTING HASS FRACTIONS AT PRIMARY 
'Z’OnE — I'N' tE'T P‘R~TN'T&'0 — AND — P-U N&H E-D — 0 U-FPUT 


IHPLICIT REAL*8 (A-H,N^Q-Z) 

NAMELIST/PRIHA/FBA, WS,RH,T 

D 0 — l-i~t — K K'K=-l-r9 ttt:: ; — :: — 

READC5?PfiTMA) 

^RlTEC^Sf PRIMA) 

NS^rT>v75-5-Sl>0- — ^ = 

02=. 231300 
AR=.0127Do 

coas-ooa-SD-o 

AIR=N2+02+AR+C02 

H2D=AIR*Ws*RH 

FUEb=TttR-*F6-A ^ 

SUm=AiR+'H20+FUEL 

N2sN2/SUM 

02 5 0 - 2 -/^ DM 

AR=AR/SUH . 

C02=Cq2/SUM 

H20SH-20-AS0M : 

FUeL=FOEL/SUM 

CH4t=FUEL 

T'F f-O'A BS (-Rh«>0-, — &Q — T-0 — 2 

WRITE CO# 100)N2,02,AR,C02f H20?CH4,T 
WPITEC7, l00)N2f 02f AR,.C02,H20»CH4, T ' 

GO — "To — lrl“l 1 

2 WRiTE(6/200)N2#O2#AR,CO2rCH4#T 
WRITEC7,2 oO)N2#02#AR»CD2,CH4,T 
0 o-P-DR-M a-T-C 0-2 =-« 5-r-'-r-AP=-’-# 0 1 1- . 5-, , - — 

UX, '-CO^^’ ,Dll,5, » ,H2Q=» / Dll,5, « ,CH4=» , D U . 5 , ' , T= » r F6 , U 

2’ rend’) 

— 2 0 0“-F0R H AT C'T'X--, J-N2st-*— f-0-2-^r#‘DT-l--» 5-f-*^ #-AR— - r D 1'1“. -#- /•# 

UX, ‘002=* ,011,5, ♦ ,CH4=‘ ,D1 U5, » ,T=' ,F6, 1 , ' &ENO' ) 

111 continoe 
end 



appendix 2 

program for CCJM.'pUTING mass FRACTIONS' AT SECONDARY 
0 fr-D-t t OT-tO N-I-N t e-T P RO &R-AH-G-I-V E:&-B OW-- PR'i-N TE B 

and punched output 


-2 ^ 5 j 6 9 1 - 0 - 

CH4 CH3 H oh H20 0 • H2 02 HCO CO 

“1-1 1-2 r3 T5 T6 

COP N mO N2 AR .N02 


-I M Pbrc-T-T— R£-A L-*r6“{TA--H'f-N>-0— ?-) ^ 

dimension COMPI (20) ,COMPDC2o) rCOMPT (20 ) , COMPE C20 ) 

C , INPUT -CONDITIONS 

N-A'M EtrI'S-T“'^-X'TNP-AH-S-rRHTF"80G>-F-&0P7-T r 

M=15 

C 

DO--l-t-H-J— K-KK--l-,-9 : 

READCSf XiNP) 

0X0=1 .OOO/FBOO 

- 0X-Pt=-1vO{>O-/PBOP- 

OXO=OXO»OXP 

COMPUTE MASS fractions OF FrESH -AIR -WITH. HUHTOlTY 

C- 1 N I T I A b-TZ E A b L— C-0 M p 0 NtN-T:S T 0 ^ -ZERO ; 

DO Hi 1 = 1, H 
UT CQMPI-(I) = 0,ODO 

COmPI C-8)'=.-23 1-300 — 

COMPI ( 1 1 ) = .0005DO 
COMPI ( 14)=.7555D0 

COMP-I ( 15-)-=.-0l-27-Da ^ ;■ 

COMPi (5)=RH*HS 

AOO = caHPI (5)+COMPI(8)+COMPI (In+C0MP-I(14)+C0MPI (15) 

COhPI (-5-) = C0HPr(-5-)-AADD - - 

COMPI (8)=C0MPI C8)/ADD 
COMPi ( Ti )=cf)Mpi ( i 1 ) / Add 

- COHPK U)=cnMPI (-14-)-/ADD 

COMPI (15)=C0MPI (I5)/AdD 



■- C OMPtH-E“frA S■S&S-OF-eoHpON6NT-S~I-^^-OI-tli■7-I•a^l— A-i-R C OM PD (-H 

DO 222 

222 COMPDf JJ3=C0MPI CJJ) *OXD 

RE-frO — !h^O'tK^rf'COH-P-T'f-I“)“,r 

PRINT 30000» (CQMPTCI) , 1=1 rH) »T 
DO 223 K=1,H 

■ - '223 C DMPT"('K~)'=C'Of+PT-EK~}'*f’a>\-P~Pi-»'ODO“l 

DO 333 K=1,M 

333 C0MP2^K)=C0HPDCK)tC0HpTCK) 

SOM-r-OT^lD’^' 

DO 444 L=UM 
444 SUm=SUH+CDMPE(L) 

■ " DO- 55 5 — 

555 COM?ecLL)=CDMPE(LL')/SUM 

WRirEc^*(-5<INPJ 

- p R j 0 000>*(-C OMPE-PI-)-, -I : — 

WRITE (7, 12345) CCGMpE Cl) , ,T 

10000 F0RHAT(70i 1 ,5,3X) 

- 12345-50RMAT-W-^,-‘-GHa=-^rD-l-5-.7 f““r&H-3^rD-b5-,7-r-^-r-H=-i-, Dl5-,-7-r-'-^ 0H=-^ 
•*SD15.7, 

1/,1X, »H2G=»,Di5,7, ' ,0=SDi 5.7, »,H2=',D15..7, »,02=’rDl5, 



2HC0=* ^Dt5.7, » ,CG=S015.7, • ,C02=» ,D15.7, =' ,D15.7, \ 

*• f /, IXr ’f^o=’ 

• • ^3-/Dl-5i7-r-'~fN2=»-,Dl-5-,7-r-*-,AR=-' rD-l5->7-f-i->— T-=« ,F6.xtf-' — • &END*- 

50000 fOrmAtC//',2X, »Ch4 = ‘ , D 1 5, 7 , 2X , ' Ch 5 = ’ , D 1 5 . 7 , 2x , ‘ H =» 
*,Dl5.7r2X, • 

— 1 OH - = D-1 5 7 , 2 X-, -J-H20--'-, D 1-5 ,-7> /-/ ,-2 Xv »-0- = -'v D 1 5 , 7 , 2X , '-H 2 - 

* = » rDl5.7,2X, 

2*02 =« /D15.7,2X, » HCO = W D15.7 ,2X, ' CO = »•, D 1 5 . 7 , //, 2X , * c 

*Q2..^,,_0i5^7r - ■■ -■ 

32X,'n =',D15.7,2X, ’NQ = ' , D i 5 . 7 , 2X , ’ N2= ' / D 1 5 . 7 , 2Xr ' AR 

* =SD15,7', //, 

4 2X,-' T--="’ >F6.-1 rV-^0 — 

iiui continue 
STOP 



